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Abstract 

Gap solitons are localized nonlinear coherent states which have 
been shown both theoretically and experimentally to propagate in pe- 
riodic structures. Although theory allows for their propagation at any 
speed V, < V < c, they have been observed in experiments at speeds 
of approximately 50% of c. It is of scientific and technological interest 
to trap gap solitons. We first introduce an explicit multiparameter 
family of periodic structures with localized defects, which support lin- 
ear defect modes. These linear defect modes are shown to persist into 
the nonlinear regime, as nonlinear defect modes. Using mathematical 
analysis and numerical simulations we then investigate the capture of 
an incident gap soliton by these defects. The mechanism of capture 
of a gap soliton is resonant transfer of its energy to nonlinear defect 
modes. We introduce a useful bifurcation diagram from which infor- 
mation on the parameter regimes of gap soliton capture, reflection 
and transmission can be obtained by simple conservation of energy 
and resonant energy transfer principles. 



* Mathematical Sciences Research, Bell Laboratories-Lucent Technologies, Murray Hill, 
NJ 07974, and Program in Applied and Computational Mathematics, Princeton University, 
Princeton, New Jersey 

t Optical Physics Research, Bell Laboratories-Lucent Technologies, Murray Hill, NJ 
07974 

■l-Mathematical Sciences Research, Bell Laboratories-Lucent Technologies, Murray Hill, 
NJ 07974 



1 



1 Introduction 



Solitons are important carriers of energy in many physical systems. The 
emergence of sohtons is understood as a consequence of the balance of dis- 
persive and nonlinear effects on the same length scale. Optical temporal soli- 
tons |l9l and more recently dispersion-managed solitons have been con- 
sidered candidates for the bits with which to transfer information over long 
distances. Recent advances in fabrication of optical fiber with microstructure 
have rendered the possibility of storing information in the form of optical gap 
solitons a natural direction for investigation. 

Gap solitons are nonlinear bound states which propagate in periodic 
structures. These have been anticipated in theoretical work |jl|, and ob- 
served in experiments ?, ^ |TB[ on sufficiently high intensity light propa- 
gation in optical fiber with a periodically varying refractive index (a uniform 
fiber grating). In contrast to bare fiber used in long distance communica- 
tions, where the formation length for solitons is on the order of kilometers, 
the formation length for gap solitons is on the order of centimeters. In theory, 
gap solitons can travel with any speed v, with < v < c, where c denotes the 
speed of light. Experiments have demonstrated the slowing of gap solitons 
to about 50% c. 

Gap solitons propagate in fibers with uniform grating structures. In this 
paper we examine gratings with localized defects to the amplitude and phase 
of the grating. We ask whether it is possible to trap moving gap solitons at 
the defect location. If so, one can envision this having important technolog- 
ical applications, e.g. optical buffers or optical memory. Through a careful 
series of numerical experiments we show how it is possible to trap gap solitons 
at a defect and elucidate the mechanism by which light energy is trapped. A 
similar question is studied in |^ using a point-particle model for the gap soli- 
ton/defect interaction. We compare our results with the conclusions drawn 
in that study. Although we refer to gap soliton capture, it is perhaps better 
called capture of gap soliton energy, for it involves the transfer of the gap 
soliton's energy to a nonlinear defect mode. 

This paper is laid out as follows. In Section ^ we derive a variable coef- 
ficient version of the nonlinear coupled mode equations from an appropriate 
one-dimensional nonlinear Maxwell model. In Section ^, we review a few 
facts about the gap soliton. We derive a family of defects which support lin- 
ear bound states (defect modes) for the coupled mode equations in Section ^. 
We then examine how these bounds states are deformed in the presence of 
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nonlinearity (nonlinear defect modes) in Section The results of this study 
are encoded in bifurcation diagrams that display the intensity as a function 
of the frequency for (a) nonlinear defect modes and (b) gap solitons. With 
the aid of these diagrams, we develop a criterion for trapping and an un- 
derstanding of its efficiency based on the notions of resonant energy transfer 
and energy conservation. Guided by this analysis, in Section ^, we perform 
a careful series of numerical experiments to show how the nonlinear bound 
states interact with the gap soliton to trap light energy. Simulations are 
carried out for the nondimensional system (|2.19|) . Dimensional experimental 
parameters are displayed in Appendix ^ Section]^ contains a brief discussion 
of the effect of nonlinear damping, a non-negligible effect in certain highly 
nonlinear materials, on soliton propagation and trapping. A summary and 
discussion of results are given in Section |^. In Appendix |^ we give physical 
parameters for silica fiber discuss nondimensionalization and tabulate the 
dimensional values of parameters corresponding to the simulations described 
in Section ^ In Appendix |B|, we describe a method for deriving defects 
supporting linear bound states with prescribed characteristics. 

2 Coupled Mode Theory in a Grating with 
Defects 

We consider propagation of light in one dimension in an optical fiber with a 
refractive index which is a spatially localized perturbation about a uniformly 
periodic index. We model the propagation of low intensity light, confined to 
a core mode of the fiber by the wave equation: 

d^[n^{z)E{z,t) ] = dl E, (2.1) 

where the refractive index is given by: 

n = n + An{^W{z) + u{z) cos {2k bz + 2<l>{z))) . (2.2) 

Here, n denotes the refractive index of the bare fiber and An, the index 
contrast, is assumed small. The functions u, $, and W model the defect and 
are assumed to vary slowly compared to the rapid sinusoidal variation of the 
refractive index. A spatially localized deviation from a uniformly periodic 
structure of period 

d = -K/ke (2.3) 
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Figure 2.1: Solid and dashed curves are two different periodic index pro- 
files with localized defects having the same "spectral characteristics" (see 
Section^. 



is obtained by taking 

u{z) ^ 1, (2.4) 
d,^{z) 0, and d^W{z) 0, as 1^1 ^ oo. (2.5) 

Low intensity light propagating in the bare fiber {An = 0) is governed 
by the spatially homogeneous linear wave equation which supports indepen- 
dently propagating, forward and backward plane wave solutions E±e^^''^~'^^\ 
where 

00 = ±-. (2.6) 

n 

The periodic structure {An ^ 0) couples these backward and forward 
components. This effect is most pronounced for wavelengths in the medium 
at or near A = 2(i or equivalently the (free space) Bragg "wavelength 

\b = 2nd. (2.7) 

For modeling propagation in the nonlinear regime we assume an instan- 
taneous nonlinear polarization 0]: 

Pnl = eoX^'^E\ (2.8) 

Combining this with ( p.2|) , the squared index of refraction with linear and 
nonlinear effects included is then 

n^{z,E^) = n'^ + nAnW{z) +2nAniy{z) cos {2kBZ + 2^ {z)) + x^^^E^. 

(2.9) 
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For high intensities the electric field evolves under the nonlinear wave equa- 
tionQ 



d^[n\z,E^) E]=c^dlE. (2.10) 

So that we can systematically obtain equations of evolution for the for- 
ward and backward carrier wave envelopes, we make explicit our assumptions 
on the medium. We assume that the variation of the refractive index is weak, 
and that the deviation from periodicity is small, i.e. there exists a small pa- 
rameter e <^ 1 such that: 

An = 0{e), 

d,W = 0{e), = 0{e), d,<l> = 0{e), and 9> = 0{e^). 

Due to the periodic structure, we expect coupling of forward and backward 
wave components. This coupling is strongest if the wavelength and period 
are chosen according to the above Bragg condition. We now make a multiple 
scales ansatz, choosing the carrier wavenumber in Bragg resonance with the 
medium: 

E = e+{z, t)e'^>'BZ+^-^Bt) + ^y-iiksz+^+^Bt) ^ (2.11) 
where the wavenumber fc^ and frequency ujb satisfy the dispersion relation 



n 



see for example [|T^]. The first two terms in (p.ll|) consist of slowly modulated 



forward and backward waves. The regime we consider is specified by the 
above assumptions on the medium and assumptions on the field amplitude, 
which we take to satisfy: 



X 



= 0{e) 



^We wo rk wit h the model (2.£)-( 2.10 ) since it yields a simple derivation of the env elop e 
equat ions ( 2.13 ). The situation is however a bit more complicated. Although ( ^.9[ )- 
( 2.10 ) incorporates the effects of photonic band dispersion, this alone is insufficient to 
arrest optical carrier shock formation on the relevant temporal and spatial scales ||l5[| . In 
fact, a valid envelope description in the absence of material dispersion would require the 
incorporation of coupling to all higher harmonics since they are in resonance. 
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The latter ensures a balance of nonlinearity and photonic band dispersion 
due to the periodic structure. We therefore anticipate that the amplitudes 
e± will be slowly varying and will satisfy: 

dte± = 0{e), d,e± = 0{e), 

The envelope functions e± in are finally determined by the constraint 

that the correction terms are of higher order in e over a time scale and length 
scale of order 

Ei/e± = 0(£). (2.12) 

The condition ( ^.12| ) requires the removal of resonant forcing terms in the 
equation for Ei. This is equivalent to the constraint that e± satisfy the 
variable coefficient nonlinear coupled mode equations: 



i-dte+ + id^e+ + Viz)e+ + K(z)e_ + r(|e+r + 2|e_r)e+ = 
c 

i-dte_ - id,e_ + V{z)e. + k{z)e+ + r(|e_p + 2|e+ne„ = 0. 
c 



(2.13) 



The coefficient functions in ( |2.13| ) are defined in terms of the parametric 
functions which characterize the index profile ( p.2|) .p| 

k{z) = ^u{z) (2.14) 
~ tt/Xt? 

V(z)^——W(z)-^'(z) (2.15) 

f=4^. (2.16) 

Our point of view is to specify grating parametric functions W{z), v{z) and 
$(z) through the constitutive law (|2.2| ). These determine the functions ki^z) 



^It is common to slightly redefine k by 

= rj ■ T:/\nv{z) / Xb 

where < 77 < 1 is defined as an overlap integral of the radial variation of the forward 
and backward modes and represents the fraction of total energy in the core of the fiber. 
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and V{z), arising in the coupled mode equations, governing the nonhnear 
propagation. Note the appearance of the combination of W and $' in ( p.l5|) . 



Therefore, spectral characteristics arising due to a DC variation (W) in the 
index can be, within this approximation, equivalently achieved through phase 
variations (<I>). The solid and dashed curves in Figure are of index profiles 
which are equivalent in this sense; see also Section W^. 



The assumptions on the variable coefficients guarantee that away from 
the defect, the system approaches the constant coefficient NLCME 

- 7T /\tI. 

V{z) ^ and k{z) ^ Roo = — • (2.17) 

We introduce typical dimensional length, Z, time, T = Zn/c, and electric 
field, S. Using these, we define nondimensional spatial and temporal vari- 
ables Z and T, and electric field E± given by: 

z = ZZ, t= T, and e±=£E±. (2.18) 

c 



Then ( p.l3|) can be expressed in nondimensional form as 



idTE+ + idzE+ + k{Z)E. + V{Z)E+ + T{\E+\^ + 2\E4^)E+ = 
idrE^ - idzE^ + k{Z)E+ + V{Z)E^ + V{\E4^ + 2\E+\^)E^ = 0, 



(2.19) 



where 

k{Z) = Zk{ZZ), V{Z) = ZV{ZZ), and T = ZS'^f. (2.20) 

Our analysis and computer simulations are carried out for the nondimen- 
sional system (|2.19|) . Conversions to dimensional form are given for impor- 



tant quantities in Appendix 0. Note that F and therefore F are positive; 
see (|2l6|). 



3 The Gap Soliton 

The linearized constant coefficient NLCME {k = Kqo, V" = 0) have a "gap" 
in their spectrum. It has no plane-wave solutions with frequencies in the 
range (— Kqo, i^oo)', see Section ^?T| . The nonlinear equations support a family 
of traveling pulses called a gap solitons 0, |^. The family is parameterized 
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by a velocity v and a detuning parameter 5 with l^;! < 1 and < (5 < vr. It 
is given by: 



E. 



sae 



2T 



A 



(sin 5) e''" sech {9-i6/2) 



where: 



—ae 
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iri 



Ho. 



2T 



(3.1) 



A (sin 5) e''" sech{9 + i6/2) 



A 



1 



= 7^00 (sin (5) (2; — vt) 



a 



1 + 

^Koo{cos6){vz - 1) ; 



a 



2 1-v^ 



e^^ + e 



s = sign(/toor) 
-id\ 



iS 



The temporal frequency of the gap soliton is KooT cos 6 which is inside the 
band gap for |f| < sin 5, although in a reference frame moving at speed v, 
the frequency is always in the band gap. We define the maximum intensity 



-^max and the total intensity /tot and give their values for the gap soliton: 



max{\E^ 

Zj 



r(3-i;2 

_ 4(1 -t;^)^ 
" r(3-t;2) ' 

and the full width at half-maximum (FWHM) by 



I tot 



1^4 



1+ l^-l 
P+ \E-\ 



KooVl - . 2 5 
sm -, 



(3.2) 



dZ 



FWHM 



2VT 



sin 5 



cosh ^ \ 1 



cos^ -. 



(3.3) 



In mathematical analysis, the square root of the total intensity /tot is often 
referred to as the norm. In Figure ^.1| , we plot the intensity as a function of 
frequency for the stationary {v = 0) gap soliton. This curve is parameterized 
by 

46 

uj{6) = Koocos5, /tot(5) = ^ foT6 e [0,7r]. 

Results on linearized stability and instability of gap solitons in different pa- 
rameter regimes are obtained in M. 
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Figure 3.1: The intensity of the stationary gap sohton as a function of its 
frequency. 

4 Defects and Linear Defect Modes 

The functions V{Z) and n{Z) define a defect in our periodic medium and we 
now seek the hnear modes associated with this defect. These are solutions 
of the hnear coupled mode equations: 



idTE+ + idzE+ + k{Z)E_ + V{Z)E+ = 
idTE_ - idzE_ + k{Z)E+ + V{Z)E_ = 0, 





or 





where 



E 





Substitution of the Ansatz 
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yields 



a,F=(j^ "(^') F, (4.5) 

\u{Z) —luj J 

u{Z) =iK{Z)e-^'^o^(0 d( (4_g) 

Solutions of the form (|4.4| ) , which are square integrable in Z are called defect 
modes. Solutions of the form ( |4.4]) which are bounded and oscillatory in Z 
are called radiation modes. The set of frequencies, lj, corresponding to defect 
modes and radiation modes is called the spectrum of ( [4 .51) . 

One can pose the question: can prescribed spectral characteristics of ( [4 .11) 
[e.g. defect modes and reflection and transmission spectra) be achieved by 
appropriate choice of u{Z) {k{Z) and V{Z))7 This question was consid- 
ered in pi| to study grating and filter design using the Gelfand-Levitan- 



Marchenko approach to inverse scattering; see also |23]. The Gelfand-Levitan 



Marchenko method, as used in ^^ , can be used to characterize gratings 



with desired spectral characteristics. 

4.1 Radiation modes for a general localized defect 

We suppose that the function u in ( |4.5| ) has the asymptotic behavior: 

u{Z) -> pe*^±, Z ^±oo 

corresponding to a spatially localized defect in the periodic structure. 

The values of u lying in the continuous spectrum are then characterized 
by the equation: 

d.F = ( %^ >""^) F. 

Seeking solutions of the form: e'^'^^v, with Q real, we find that there are 
nontrivial solutions v provided: 

= + (4.7) 

Therefore, the continuous spectrum consists of the real axis minus a gap 
(photonic band gap): —p < u < +p; see Figure 
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4.2 Dark soliton defect gratings 



It is interesting to note that the system ( [4 .51) , characterizing the modes of 
the modulated periodic structure, is the Zakharov-Shabat eigenvalue problem 
associated, via the inverse scattering transform, with the defocusing nonlinear 
Schrodinger equation p5|: 



idrU - d^zu + 2|m|V (4.8) 

If u{Z, t) satisfies (|4^ ) then for each r, (a) the spectrum of (|4.5|) is indepen- 
dent of r and (b) as r varies, the eigenfunctions and radiation modes of ( [4 .51 ) 
evolve in a trivial (linear) and explicitly computable manner. That is, equa- 
tion ( |4.8| ) defines an isospectral deformation of the eigenvalue problem ( [4.5| ) 
and therefore, by ([4.6|) , provides a rich class of potentials, u{Z), and therefore 
modulated gratings {k{Z),V{Z)) with the same spectral characteristics. 

From the inverse scattering theory of ( [4.8| ), we learn that the bound states 
and continuum radiation modes associated with ( ^^ ) can be "mapped" , re- 
spectively, to the dark solitons and radiation modes of ( [4.8|) . In this section 
we explore those grating profiles which correspond to the simplest dark soli- 
ton solutions of ( |4.8| ). 

Since the eigenvalue problem for the pair {F{Z),uj) is self-adjoint, uj varies 
over the real numbers. Furthermore, the set of all u satisfying |ti;| > p is 
continuous spectrum. There are no eigenvalues embedded in the continuous 
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spectrum. Therefore, if the eigenvalue problem ( |4.5| ) has eigenvalues they 
must occur in the gap \uj\ < p. Of interest is the following 
Inverse Problem: Given N numbers a;i, . . . ,ujni satisfying \ujj\ < p, find 
potentials, u{Z) for which these are the eigenvalues of the eigenvalue prob- 
lem (p|). 

This inverse problem has many solutions. A remarkable class of solutions 
are those, for which the reflection coefficient associated with ( [4.5|) is zero. 
These are the dark N-soliton solutions of defocusing NLS. 
Dark solitons (N=l): Let 7^ be arbitrary and Ic^l < p. Let 

p = \uj + ik\. 

Define 

u{Z) = e''^{yj - ik tanh {kZ) ) (4.9) 

where is left unspecified to this point. 

It can be verified easily that A is an eigenvalue of ([4.5|) with corresponding 
eigenfunction: 

F{Z) = (^[^^)sech{kZ), (4.10) 



Therefore, E± are given by: 



i arctan ^iiSHMfcZ) 



, , , , i J ktan.h.(kZ) 

El V ze~^^e~ 2 ^^^^^^ — — 



6-*'"*sech(A;Z). (4.11) 



Specification of the grating with prescribed defect energy: 

A grating is specified by the functions V{Z) and k(Z) (see ( p. 151) and ( p. 14]) 
and (|2.20|) ). Using (4^) we obtain a relation between a family of "dark soli- 
tons", u{Z), and the functions V and k: 

iK{Z)e-^'^o^(')d' = {uj - ikiaxi\i{kZ)) . (4.12) 

Choosing = 7r/2, this yields 

k{Z) = [iu^ + tanh^ (kZ)] ^ , (4.13a) 
V{Z) = l-P Lu [uj^ + k^ tanh^ (kZ)] sech^ {kZ) (4.13b) 
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and sets the term 

ie-'^ = 1 



in equations ( |4.10| ) and ([4.11|) . Note that the hmit uj of the defect 



definition (|4.13| ) is singular: 

k{Z) = \ktanh{kZ)\, (4.14a) 
V{Z) = ±^6iZ) (4.14b) 

where S{Z) denotes the Dirac delta function and the sign on V{Z) depends 
on the direction on which the limit ^ is taken. Instead, taking V{Z) = 0, 
uj = 0, e**^ = ±1 in equation ( [4.12 ) yields the continuous limit 



k{Z) = -ktanhkZ (4.15) 
F{Z) = sech kZ (4.16) 

Recall that in our coupled mode system, the medium is characterized by 
three functions: z/, W, and $. From our solution we see by ( |2.14|) that (|4.13a|) 



(or (|4.15|) ) uniquely determines u. However, from (|2.15|) we see that ( gl3bD 



specifies only a linear combination of W and giving one some freedom 
in how to design a medium with the desired spectral characteristics. In Fig- 
ure two gratings with identical z/, k, and V are displayed. The solid curve 
corresponds to the choice $ = (no phase shift in the the refractive index), 
and the dashed curve corresponds to the choice W{Z) = (no modulation 
to the DC component of the refractive index). 

4.3 More general defect gratings 

In this section we consider a class of defect gratings which generalizes those 
studied in the previous section: 



k(Z) = '^cj2 + n2ptanr (A;Z); (4.17a) 

ViZ) = ^^k^^ech'(f) 
^ ^ 2{uj^ + n^kHanh'{kZ)) ^ ' 

We shall use this class of defects extensively in numerical simulations. This 
family of defects can be obtained by specifying an exponentially localized 
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Figure 4.2: The three eigenmodes for the defect ( [4.17[ ) with (uj,k,n 



(—1,2,2). In each plot is in the left column and E_ in the right. Solid 
and dashed lines correspond to real and imaginary parts, respectively. 



eigenfunction of ( |4.5| ) and then deriving a potential for which this function 
is a bound state. The calculation is presented in Appendix 0. The general- 
ization of the eigenmode ( [4.11| ) 



IS 



^ arctan ^^^^ C^-^) 



i . nk tanh (kZ) 

E^J V e~ 2 ^ — 



e-*'^*sech"(fcZ). (4.18) 



When n > 1 in ( [4. 171) , then the linearized equations ( [4.1|) have multiple bound 
states. For example when {uj,k,n) = (—1,2,2), there are three eigenvalues 



Uq = —1 and u±i = ±v 13. The eigenfunctions are shown in Figure |472 . 

The photonic band gap width is determined by the asymptotic behavior 
of K in (|4.17a| ) and is given by 2k,oo, where 

Koc = lim k{Z) = Vuj'^ + n'^k^. (4.19) 

In our computer simulations (see Section P) we will set k'^ = k"^ = C'^/n'^, 
thereby fixing the gap width equal to + and vary the defect width 

by varying the parameter n. 

The defect k{Z) varies between the values = Ic^l and Kqo = V<^^ + n'^k'^. 
Therefore the "depth" A* of the defect is given by 

A, = + n^k^ - \iu\. (4.20) 
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We define the defect widtli (FWHM) to be given by twice tlie value of Z 
for which = + Kqo), yielding 



2 / J2|c<;|V^j2qr7?F + n2p-2cu2 

FWHM = -tanh"^ I (4.21) 

h \ Ink I V ; 

in the nondimensional setting. Dimensional values for the defect depth and 
width are provided in Appendix 0. 



5 Nonlinear defect modes 

In this section we show that the linear defect modes of Section upon 
inclusion of nonlinear terms, deform into nonlinear defect modes in a sense to 
made precise below. We begin by observing that the dimensionless nonlinear 
coupled mode equations ( |2.19|) can be written in the vector form: 

[% {dr + (yj,dz) + + V[Z))E + VN{E, E*)E = 0, (5.1) 

where E is the two- vector with components E±, ai and are displayed 
in ( |4.3| ), and the N{E,E*) defines the nonlinear term: 

iv(E,£.,^(i^-i^;2|^-i= (5.2) 

Bifurcation of nonlinear defect modes We assume that all eigenvalues 
of the linearized problem ( [4.1D are simple, which has been found numerically 
for the family of defects investigated. Let Eq = e~'''^°^ So{Z) denote a linear 
defect mode. That is, So{Z) is a spatially localized solution of the equation 

{uJo + icT3dz + cT,K{Z) + V{Z))E = 0. (5.3) 

We seek to construct nonlinear bound states of ( p.lj ) of the form: 

E{Z,T) = e-''^^ S{Z), (5.4) 

where 

£iZ) = a{EoiZ) +\a\^E,{Z) + Oi\a\^)) 



(5.5) 
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and a is a small parameter. Since for any cj, S{Z) = is a solution of ( p.3| ), 
solutions of the form ( ^.Sj) are said to bifurcate from the trivial solution at 

UJ = UJo- 

Substitution of ( ^.51 ) into (|5.1| ) yields a hierarchy of inhomogeneous linear 
equations beginning with: 

: Co Eo = (5.6a) 
0{\a\') -.CqE, = -UJ, E, - T N{Eo,E*)Eo (5.6b) 

where 



Co = ujq + ia^dz + c^in + V. (5.7) 



is a linear self-adjoint operator. The first equation in (|5.6|) holds if Eq = Sq, 
any linear defect mode, So of frequency uq. The eigenvalue uo is necessarily 
of multiplicity one. The second equation in ( |5.6|) is solvable for a localized 
correction term, Si if and only if the right hand side of the equation is 
orthogonal to the null space (zero energy subspace) of Cq. Imposing this 
orthogonality condition yields the following equation which determines the 
value of Ui. 



{So \ujiSo + T N{So,S;)So) = 0. (5.8) 



We obtain from (|5.8| 



^ (£o|iV(£o, £„•)£„) 
= (£„!£„> <^-'" 

It follows that the nonlinear defect mode bifurcating from the linear defect 
mode of frequency oJo is: 

E{Z,T) = e-*"^« {So{Z) + \a\^Si{Z) + 0{\a\')) (5.10) 
u; = u;o+|tt|V + C'(|a|^) (5.11) 

We shall refer to the nonlinear defect mode bifurcating from o; = cuq as an 
cuo-nonhnear defect mode. Since F > (see Section the bifurcating states 
have frequencies below Uo- 

A rigorous proof of the existence of bifurcating nonlinear defect modes 
and the validity of this expansion can be given in a manner analogous to 



that carried out in the context of the nonlinear Schrodinger equation |20 . 
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Experiment 1 .2 


- 




^S,,^ Experiment 1.1 











Frequency 

Figure 5.1: Intensity vs. frequency for the gap soliton (bold) and a nonlinear 
defect mode with parameters {uJo,k) = (—1,4). 



Numerically we are able to find nonlinear defect modes for values of a much 
larger than zero. A plot of the intensity vs. frequency for one such family 
of defect modes is shown in Figure |5.1|. For the spatially homogeneous case, 



gap solitons are seen to bifurcate from the zero state at the right endpoint of 
the continuous spectrum, u = = \/l7. For the given defect, a branch of 
nonlinear defect modes bifurcates from the zero state at u = ujq = —1. The 
nonlinear defect mode and its frequency become difficult to compute as the 
frequency u of the nonlinear mode approaches an endpoint of the continuous 
spectrum, because the exponential decay rate decreases and larger spatial 
intervals must be used in order to compute the exponential tail. 



6 Computer simulations of a gap soliton in- 
cident on a defect 

In Section | we discussed gap solitons, the fundamental nonlinear bound 
state of propagation in a uniform periodic structure. In Section ^ we then 
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considered the linear modes of a periodic structure with a defect and, in 
Section ^, the nonhnear defect modes which bifurcate from these hnear defect 
modes. In this section we study by computer simulation the dynamics of a 
gap soliton incident on a defect. 

In Section ^, we derived a nondimensional form of the coupled mode equa- 
tions with nondimensional parameters n{Z), V{Z), and F given by ( p. 20 ). 
By an additional rescaling of the nondimensional variables Z, T, and E±, 
we may redefine the equations so that F = 1 and Koo takes the form given 
in equation (^4.19|) . In other words, Kqo and F can be set arbitrarily and the 
results can always be mapped to a physical system. 

Simulations indicate complex interactions between the incident gap soli- 
ton (Section 1^) and the modes of defect (Sections ^ and |). An understanding 
of the dynamics and the potential for trapping requires an understanding of 
the energy exchange between the gap soliton mode and nonlinear defect modes. 
Our numerical simulations give strong support to the following hypothesis 
suggested by the notions of resonant energy transfer and energy conservation. 

Hypothesis. Consider a gap soliton incident on a defect with sufficiently 
low incident velocity. The gap soliton will transfer its energy to a nonlinear 
defect mode, and thereby be trapped, if there exists a nonlinear defect mode 
of the same frequency and lower total intensity (L"^ norm). Otherwise, the 
gap soliton energy will be reflected and/ or transmitted. 

We now describe our numerical experiments. The defects described in 
Section ^ are members of a 3-parameter (cuq, fc, n) family and the gap soli- 
tons are described by two parameters (f, 5). While it is not possible to 
investigate all of the soliton-defect interactions in this 5-dimensional space, 
we have performed a large number of simulations, and were able to draw some 
general conclusions. We concentrate on gap solitons of comparable width and 
amplitude to the defects, so that linear and nonlinear interactions are likely 
to be strong and balanced. Physical experiments have so far produced pulses 
with relatively large values for v and small values for 5, so we make some 
attempt to trap gap solitons in this parameter regime. Dimensional equiva- 
lents for most of the nondimensional experiments are given in Appendix 
Intensities range from 130 to 1800 GW/cm^, pulse widths are between 1.3 
and 4.4 mm, and defect widths are between 1.6 and 4.7 mm (FWHM). Note 
that in this section, the frequency of the linear defect mode is given by cuq, 
while uj = Koo cos 5 is the frequency of a stationary gap soliton. 
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Figure 6.1: (Experiment 1.1) The initial value of l-E+P + l-E-p (solid line), 
which gives the approximate strength of the nonlinear forcing, and of the 
defect k{Z) — (dashed), which gives the forcing due to the defect. 



6.1 Experiment 1: Gap soliton incident on a dark soli- 
ton defect gratings 

We first consider the simplest case of pulses interacting with the dark soliton 
defect gratings defined in ( [4.13|) . We consider a defect with k = A and ujq = 



— 1. This defect supports a single nonlinear bound state. The key insight 
into predicting whether the gap soliton interacts strongly with the defect is 
found by examining Figure ^.1| . A gap soliton with frequency u = cos 5 
will interact most strongly with the defect if a nonlinear defect mode exists 
with the same frequency and equal or less total intensity. If it does, then it 
is possible for the gap soliton to resonantly transfer its energy to the defect 
mode. The relevant mechanism seems not to be the slowing of the soliton, 
rather this transfer of energy. Note from the figure that for uj < —1 {5 > 1.82) 
such modes exist, and for larger uj they do not. Of course, the gap soliton 
frequency-intensity curve given in the figure applies only in the case f = 0, 
but it is useful for making predictions for small v. 

Experiment 1.1 (Reflection/Transmission) With the defect parameters set, 
we then investigate a two parameter family of gap solitons indexed by the 
velocity v and the detuning 5. Our first experiment is for detuning 6 = 0.9. 
Via inspection, the pulse is of comparable depth and width to the defect, so 
seems an ideal candidate for capture, see Figure lO] . However, the central 
frequency of this gap soliton at small velocities is Koo cos 5 ~ 2.56 > ujq = — 1, 
so we do not expect the defect mode to be strongly forced by the gap soliton. 
Indeed, although we observe a slowing, and therefore delay, of the gap 
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Figure 6.2: (Experiment 1.1) At left, a gap soliton with v = .2565 is reflected 
by a defect ai Z = 0. At right, a slightly faster gap soliton with v = .257 is 
transmitted. 



soliton we do not find significant excitation of the defect mode or trapping. 
We find that below a critical velocity v ~ 0.257, all gap solitons are reflected, 
and above this speed they are transmitted. The closer the incoming pulse 
comes to this incoming velocity, the longer it remains in the neighborhood 
of the defect before being ejected and the velocity of the outgoing pulse is 
approximately that of the incoming pulse. 

In Figure |6.2| , we show the evolution of two gap solitons incident on the 
defect, both very close to the critical velocity, showing clearly the effects 
discussed above. 

Interestingly, the gap soliton slows down when it nears the defect (Fig- 
ure |6]^ for times between 20 and 30). This is somewhat unexpected; as the 
defect supports a bound state, one intuitively expects the soliton center of 
mass to move as a "classical particle in a potential well" . Instead the soliton 
behaves more like a classical particle encountering a potential barrier. Brod- 
erick and de Sterke conjecture that if a defect supports bound states, then 
a particle approaching it should "see" a potential well, and if it supports no 
bound state, then an approaching gap soliton should see a barrier Our 
numerical simulations and theory, based on energy conservation and resonant 
energy transfer illustrate that the situation is more complex. Indeed, both 
the "potential well like" and "potential barrier like" behavior are possible for 
a defect that supports bound states. 
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Figure 6.3: (Experiment 1.2) With 5 = 2, a. defect mode of significant inten- 
sity remains behind after the sohton passes through. (Due to the absorbing 
boundary conditions used in the simulations, the gap sohton dissipates as it 
approaches the edge of the computational domain) 



The reflection of these gap solitons is well-explained by Figure 5.1. To 



test our hypothesis, we post-process the numerical experiment as follows. At 
each time step, we compute the projection of the solution onto the linear 
defect mode of this defect, and find that when the gap soliton is directly over 
the defect, the projection onto the bound state accounts for only 6% of the 
total norm of the solution, and after the soliton escapes, the projection 
accounts for less than 0.2% of the solution. As there is no resonant exchange 
of energy between soliton and defect mode, the soliton escapes. 
Experiment 1.2 (Trapping for larger intensities) Figure ^]l| suggests that 
gap solitons with larger values of the detuning S, may interact more strongly 
with the defect. We therefore run the experiment again with 6 = 2 (fre- 
quency u = —1.72) below ujq, and v = 0.2 the results of which are shown 
in Figure |6.3| . When the soliton encounters the defect, it seems to split into 
three parts: a transmitted soliton, a trapped mode, and radiation. The mode 
that remains at the defect has only about 16% of the total intensity of the 
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incoming gap soliton. Remarkably, at the end of the computation, the cap- 
tured state's frequency is approximately u = —1.7 and the solution's total 
intensity is such that the trapped energy is described by a point very close 
to the nonlinear bound state curve of Figure This supports the first 
part of our hypothesis. For small values of v, the amplitude and frequency 
of the trapped state do not seem to depend on v. Above a certain larger 
velocity, significantly less energy is trapped by the defect, suggesting that 
the hypothesis needs refinement for large velocities. 

Experiment 1.3 (Refining these results) Clearly, this computation shows 
that we can use a defect to trap a significant portion of the electromagnetic 
energy in a gap soliton. However to this point, the gap solitons we have 
captured have had very high intensities, and further, only a small amount of 
the solitons energy is trapped by the defect. It is of interest to trap lower 
intensity pulses, and it would be preferable if a larger fraction of the soli- 
ton's energy were trapped by the defect. The nonlinear defect modes always 
bifurcate to the left from the linear defect mode frequency for increasing 
intensity. Consider the defect defined by ( [4.13[ ) with /c = 4 as before, but 



now letting uq = 1. This leaves k{Z) unchanged, while changing the sign 
of V. This moves the base of the nonlinear bound state curve to = +1, 
so that the intensity-frequency curves for the nonlinear defect mode and the 
gap solitons are significantly closer together; see Figure |6.4]. If we examine 



the interaction of 5 = 7r/2 (cj = 0) gap solitons with each of these defects, 
the bifurcation diagrams anticipate that the Uq = 1 defect will capture a lot 
of energy from the pulse, while the ujq = —1 defect will reflect or transmit 
the pulse, depending on its incoming velocity. Numerical experiments show 
this to be the case. 

We can further improve trapping using the dark-soliton family of defects 
by increasing the ratio ujq/koo- The gap edge is at k^o = \/uJq + /c^, while the 
defect mode curve starts at Uq and goes left with increasing intensity. The 
problem with this approach is that as that as Uq Koo, the the width of the 
gap increases without bound, while its depth goes to zero. 

6.2 Experiment 2: Gap soliton incident on grating de- 
fects supporting multiple bound states 



We now show how to use the generalized dark soliton defects of Section 4.3 



to more efficiently capture solitons. As pointed out in the paragraph follow- 
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Figure 6.4: Intensity vs. frequency for the gap soliton (bold) and a nonlinear 
mode for defect with parameters (cuq, k) = (1,4). Although the defect is the 
same as that used in Figure |5.1| excepting a sign in the definition of V{Z), 
the defect mode curve is further to the right and closer to the gap soliton 
curve, predicting greatly improved trapping. 
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ing ( ^4.19| ), we can use defects given by ( [4.17| ) to fix the spectral gap and study 
the interaction of gap sohtons with defects of different widths. We next study 
gap sohtons incident on a grating of this form with {LjQ,k,n) = (—1,2,2). 
This defect is twice the width of the dark sohton defect grating of Section |6.1| , 
but has the same hmiting profile far from the defect region. From, Ap- 
pendix 0, the defect supports three hnear bound states, with ground state 
frequency luq = —1 and excited states u;±i = ±vT3. Branches of nonlinear 
defect modes bifurcate from each of these linear modes. Figure |6.5| is the 
analog of Figure ^]T] for this defect. To the left of the indicated frequency 
cj*, the cu+i-nonlinear defect mode has greater intensity than the gap soliton. 
The frequencies u^i, Uq, tu*, and cj+i divide the band gap into 5 regions. 
In the regions — Kqo < uj < uq and < u < u^i we expect by the same 
mechanisms as in Experiment 1.2 to find trapping of energy, while in the 
regions Uq < uj < and u^i < oj < Koo, we do not expect trapping. In the 
segment —Koo < uj < uj_i, we expect complex behavior because two trapped 
nonlinear modes coexist. We expect the most efficient capture for solitons 
with frequency slightly greater than tu*, for which a nonlinear bound state 
exists of slightly lower total intensity than the incoming pulse. 
Experiment 2.1: (The trapping region cu* < < cu+i) 
(2.1a) We first examine gap solitons with 6 = 0.9, {u ~ 2.6) which lie just 
to the right of cj* in Figure |6.5| . Trapping here is relatively efficient, because 
a soliton can transfer almost all its amplitude to the nonlinear defect mode 
of the same frequency and slightly lower intensity. We find trapping for gap 
solitons slower than a critical velocity of about Vc = 0.102. In Figure |6.6| , 
we show the evolution of a gap soliton, initially propagating to the right, 
which gets trapped. In Figure |677| , we show the position vs. time plot for a 
gap soliton that gets trapped and one that escapes. In both cases, the gap 
soliton speeds up on reaching the defect, consistent with the defect acting as 
a potential well. Figure |6]^, shows the output soliton velocity as a function 
of the soliton input velocity in a region near the critical velocity. The figure 
indicates a sharp transition at a critical velocity from gap solitons which are 
trapped to gap solitons which propagate through. 

Although the NLCME system ( pA9D conserves /tot; radiation may carry 
some energy away from the defect. The computations are performed on a 
finite domain with absorbing boundary conditions, so that radiation losses 
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Figure 6.5: Intensity vs. frequency for the gap soliton (bold) and the three 
nonhnear mode for defect with parameters {ujQ,k,n) = (—1,2,2). Trapping 
is possible for frequencies on the thickened section of the gap soliton curve. 



Figure 6.6: (Experiment 2.1a) A typical picture of the capture of a gap soliton 
by a defect centered at Z — 0. 
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Figure 6.7: (Experiment 2.1a) The position vs. time of an escaping and a 
captured gap soliton. Note that the instantaneous velocity increases when the 
gap sohton is in the defect region [Z near zero). 6 = 0.9, v near Vc ~ 0.103, 



the potential as described in subsection 3.2 
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Figure 6.8: (Experiment 2.1a) Initial velocity (fj) vs. final velocity {vf) of 
gap soliton, parameters as in with variable v. 
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Figure 6.9: Decay of the local norm (normalized) for Experiment 2.1a. 



can be measured by monitoring the local norm, i.e. 

Local norm = (^j \E+\^ + \E_\^dZ^ ' , 

where D is a bounded region containing the defect; see Figure p.9| . By time 
t = 120, the energy has been transferred from the soliton to the nonlinear 
defect mode, but the system continues losing energy to radiation at a constant 
rate for the length of the simulation. 

(2.1b)For the slightly a smaller value of 5 = 0.6, the distance between the 
gap soliton curve and the nearby defect-mode curve in Figure ^3] is increased. 
Some of the energy is deposited in a defect mode, while the remaining energy 
appears to propagate as a diminished gap soliton plus small radiation; see 



Figure 3.10 



There is not space to report in detail on all of the behaviors found in 
the numerical simulations for this defect. As expected, gap solitons with 
frequency uj < ujq are trapped in a similar manner to those in Experiment 1.2. 
When uj < uj^i, then all three nonlinear modes are excited by the gap soliton. 
More unexpectedly, in the region c^o < < t^*, the gap soliton, while never 
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Figure 6.10: (Experiment 2.1b) Partial capture; an incident gap soliton re- 
sults in part of its energy captured in the defect and part transmitted as a 
lower energy gap soliton 
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Figure 6.11: (Experiment 2.2a) Capture of a gap soliton with 6 = 0.45 by a 
wide well. 




Figure 6.12: Local norm of the solution (bold), the projection onto the 
eigenmode (solid) and the cj+i mode (dashed) for Experiment 2.2a (left) 
and 2.2b (right). 



captured, is never reflected either. In this frequency range, for every initial 
velocity as low as f = .0006, the soliton is transmitted after slowing down 
slightly when encountering the defect. 
Experiment 2.2: (Wider defects) 

(2.2a): By widening the defect, we may place more eigenvalues closer to the 
edges of the band gap which might then be used to trap gap solitons with 
even smaller 6. Using a defect with parameters {uJo,k,n) = (—1,1.6,2.5) 
(keeping Kqo = VTf as in previous sections), which has 5 eigenvalues c^o = 
— 1, uj±i = ±-\/281/5, uj±2 ± V409/5, we captured a soliton with 5 = 0.45, 
although with the velocity significantly reduced to about v = 0.025. We 
found that the defect with k = n = 2 described in Experiment 2.1 reflects 
this gap soliton, as its central frequency is to the right of the defect mode 
curve in Figure |6.5|. The dynamics of the soliton captured by the present 



defect is shown in Figure |6.11| . At about t = 400 the trapped mode begins to 
lose intensity to radiation, as is more clearly seen in Figure |6.12| a; note the 
decay of the local norm beginning around t = 400. 
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Also shown in Figure 6.1!j are the numerical projection onto the linear 
eigenfunctions belonging to uj^2 and uj^i. At capture time (t ^ 60), the 
solution is dominated by the ui+2 mode. At longer times, however, this 
trapped mode is not persistent and the energy is transferred from a;+2 to 
tu+i, with a lot of energy lost to radiation and very little in the other three 
eigenmodes. 

(2.2b): The effect of defect width can be further explored. In the next 
experiment, we widen the defect further, by choosing parameters k = 4/3, 
n = 3 in (|4.17a|) , and choose the same gap soliton parameters as in the above 
paragraph. The defect is slightly wider than in Experiment 2.2a, though it 
still supports five linear bound states, with eigenvalues slightly smaller than 
in the previous paragraph: o^o = —1, uj±i ± \/89/3, uj±2 = ±v^l37/3. In this 
case, the trapping is significantly less effective. A much smaller bound state 
is trapped. As in the previous experiment, most of the energy is localized in 
the modes belonging to tu+i and a;+2, although, as seen in Figure |6.12| , the 
tu+i-mode begins growing sooner. 



6.3 Experiment 3: Defect arrays 



Figure ^]8| shows that gap solitons which are not trapped by a defect may 
be severely slowed. This suggests that an array of defects of the type dis- 
cussed above can be used to successively slow and then trap a gap soliton. 
Using a pair of defects, we have been able to trap a gap soliton whose initial 
velocity was 50% higher than the critical velocity found in Experiment 2.1. 
We construct defect arrays k.2{Z) and V2{Z) by forming k.{Z) and V{Z) 
of subsection |6.2| and then letting k,2 = i^{Z — Zi) + k,{Z — Z2) — Koo and 
V2iZ) = V{Z - Zi) + V{Z - Z2). Such gratings with Zi = -3 and Zi = 3 
are shown in Figure |6.13| . In Figure |6.14| we show the position vs. time for 



gap soliton with S = .9, v = .15, which is slowed by the first well and then 
captured by the second. 



6.4 Experiment 4: Side barriers 

In the previous sections, we captured light by transferring energy from a 
gap soliton to a nonlinear defect mode. In this section, by contrast, we trap 
a moving soliton between two obstacles. In Experiment 1.1, solitons with 
frequency to the right of the ground state frequency were slowed but not 
trapped upon encountering the defect. We modify the defect by adding a 
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Figure 6.13: (Experiment 3) k{Z) and V{Z) for an array of two defects. 
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Figure 6.14: (Experiment 3) Position vs. time for a gap soliton incident on 
an array of defects. Defect positions given by dashed lines. 
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Figure 6.15: (Experiment 4) The modified defect described in 



6.4 



bump or "potential barrier" away from the main defect; see Figure |6.15 . 
This configuration of defects traps the gap sohton in a novel way. Instead of 
a bound state forming near the minimum of k{Z), the gap soliton bounces 
back and forth between the old "well" and the new "bump" that has been 
added. Further, it captures a pulse with an incident speed of f = .3, about 3 
times the critical velocity for the generalized dark soliton grating described in 
Experiment 2.1. In addition, the rate of energy loss for gap solitons captured 
by this defect is significantly reduced. Figure |6.16| shows the rate of energy 
loss for gap solitons with 6 = .9 and velocities v = .2625 and v = .3. It 
shows that, although it can capture gap solitons moving this fast, as the 
speed increases, the efficiency of the capture decreases. 



7 Nonlinear damping effects 

As we have shown in the preceding section, the gap solitons for which we 
have been able to find interesting capture behavior have all had large values 
of 6. (So large, in fact, that perhaps silica fibers would not be able to support 
pulses of that intensity.) One potential way to make use of our theoretical 
solutions would be to use fibers with larger nonlinear refractive index n2. 
We see by equation ( |A.9| ) that the intensity of a gap soliton, for fixed v and 
6, is inversely proportional to n2. Chalcogenide fibers, for which 77,2 is as 
much as 500 times larger than in silica fibers, are a promising material in 
which one could potentially observe the phenomena discussed above at lower 
intensities. Unfortunately, in chalcogenide fibers, the imaginary part of n2 
is also significantly larger, corresponding to nonlinear damping arising from 
multiphoton absorption. By chooseing a chalcogenide glass composition that 
minimizes both two photon and three photon absorption, one can achieve 
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Figure 6.16: (Experiment 4) Local norm as a function of time for gap 
soliton with v — .2625 (solid) and v — .3 (dashed). 
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an n2 nearly 500 times silica while suffereing a multiphoton loss of of a few 
percent at intensitiess requeired for a nonlinear phase shift of tt. 

In terms of the coupled mode equations, a complex cubic refractive index 
gives rise to a complex coefficient T. Due to the symmetries of the NLCME 
and the gap soliton, the magnitude of F is unimportant, as gap solitons have 
intensity that scales as 1 /F. What will be important is the ratio Fj/F^. In this 
case we are more interested in simply simulating the propagation of pulses 
which at t = correspond to gap solitons. The strength of the damping is 
proportional to /max, or sin^ |. Therefore as we decrease 6 the effect of the 
nonlinear damping should be decreased. However, if decreasing S requires a 
decrease in the gap soliton velocity v for trapping, so that pulses will have 
more time to decay as they propagate before reaching the defect. We ran 
one set of experiments with Fj/Fr = 0.1, 0.01, and 0.001, with v = 0.2 and 
6 = 0.9, and with the generalized dark soliton grating with ujq = —1, k = 2, 
and n — 2. Without nonlinear damping, the defect will not capture these 
solitons. The soliton was initialized 5 units to the left of the defect center. 
With the damping ratio .1, the gap soliton is mainly damped before it even 
reaches the defect. With ratio .01, the pulse loses just enough energy that 
much of it is captured upon reaching the defect. With the ratio .001, the gap 
soliton propagates through the defect untrapped. 

8 Summary and discussion 

Gap solitons are localized nonlinear bound states which propagate in periodic 
structures. We have investigated by analytical and numerical methods the 
possibihty of capture of gap solitons by the introduction of appropriately 
designed defects, spatially localized deviations from exact periodicity. 

We first displayed interesting classes of defects which support trapped 
defect modes of the linear coupled mode equations. We then showed that 
these linear defect modes, deform into nonlinear defect modes of the nonlin- 
ear coupled mode equations. Bifurcation diagrams of total intensity (/tot) 
versus frequency suggest the following hypothesis: for sufficiently low veloci- 
ties, a gap soliton incident on a defect will transfer its energy to a nonlinear 
defect mode localized at the defect provided there is one of the same fre- 
quency (resonance) and lesser total intensity (energetically accessible). This 
hypothesis is supported by an extensive series of numerical investigations. 
An understanding of the dependence of the critical velocity requires further 



35 



investigation. 

We have studied the interaction of gap sohtons with defects which support 
one or multiple linear (and therefore nonlinear) defect modes. In Section |6]^ 
we show how secondary defect modes can be used to trap energy from lower- 
amplitude gap solitons, which correspond more closely to the regime accessed 
thus far in physical experiments. However, in contrast to trapping by a 
defect mode with a single defect, when gap soliton energy is transferred to 
the modes of a multimode defect, the dynamics of the localized energy is 
quite complicated. 

We believe that a finite-dimensional model incorporating both soliton and 
defect mode degrees of freedom could be very useful in understanding the 
capture problem. In these models the soliton is modeled by several param- 
eters {e.g., position, width, phase) and the defect mode by its "amplitude". 
A system of ordinary differential equations approximating the dynamics of 
a soliton interacting with a defect can be obtained from an effective La- 
grangian, which is a function of these collective variables. Broderick and 
de Sterke ^ have studied such a model which does not take into account 
degrees of freedom available in the defect modes. Their model displays some 
of the observed behaviors but, as we have seen, the mechanism of resonant 
energy transfer must be included to provide a full explanation. 

Similar finite-dimensional models have been studied in [|r^, [TT], |T2|] for 
the sine-Gordon, 0^ and nonlinear Schrodinger equations. Comparison of 
models for the trapping of sine-Gordon kinks with and without the defect- 
mode degrees of freedom |ll| demonstrates the necessity of allowing the 
additional modes of oscillation. We have applied tools of dynamical systems 
theory to similar reduced models for soliton like structures of the sine-Gordon 
equeation and NLS [|13| . These studies give insight into the nature of the 
set of states incident on the defect resulting in transmission without capture, 
capture, and capture for all time. Closer qualitative agreement with the full 
dynamics is obtained by inclusion of a damping term, reflecting the coupling 
to radiation modes. 



A Calculation of gap soliton characteristics 



In this appendix, we derive dimensional values of Koo and F in equations (|2.14|) 



and ( |2.16| ). We use these quantities to estimate the physical parameter values 
corresponding to the simulations performed in Section ^ 
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A.l Dimensional values of koo and F 

It is common to work with the cubic refractive index [0, p. 40, 582], n2 (also 
denoted nl): 

n2 = 3x^^V4eocn', (A.l) 
quoted in units of m^/W. This gives 

~ 47reocnn2 



A, 



(A.2) 



We hst the parameters needed, as obtained from Eggleton, et. al. and 
from standard sources: 



\b = 


1053 nm; 


n = 


1.45; 


An = 


3 X 10"^; 


n2 = 


2.3 X 10-2° mVW; 


c = 


2.98 X 10®m/s; 


eo = 


8.85 X 10-^2 C/Nm^ 



(A.3) 



Then the wavenumber in the medium is given by 



71/ 

kB = 27T— = 8.7 X 10^ m-\ (A.4) 
As 



and Koo (see (|2.14|) ) and T by 

«:oo = 900m-\ and (A. 5) 

f = 1.06 X 10-1^ (A.6) 
N^m 

A.2 Converting from nondimensional to dimensional 
form 

Inverting the nondimensionalization relations ( |2.20 ), the dimensional length, 
time and electric field scales can be obtained from the following relationship 
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Koo TfAn ' nAn c 

£^ = 2^ = T^^, (A,7b) 



where we have also used (^.141) and (|2.16|). 



A. 3 Dimensional Experimental Parameters 

In terms of the solution to the NLCME, the (scalar, dimensional) electric 
field is given by 

where cc represents the complex conjugate. The mean amplitude of the 
electric field is thus given by 

\E\^ = 2 (|e+|2 + |e_p) = 28^ (l^+p + \E.\^) , 

neglecting phase and cross terms. The maximum intensity |0| is 

I = -eocnmaxli?!^ = eocnS^Imayc- (A. 8) 

Note that this is a scaling of the nondimensional quantity /max defined 
in (^). Combining this with and (PTTED , 



722(3-1;^) 2 



Scaling ( p.3|) by (|A.7a|) and using that Kqo = VLJ^~l^ri^ gives 



Dimensional FWHM = — ■ cosh ^-1/1 + cos^ -. 

TrAn sino V 2 

The temporal width is just 



Dimensional FWHMtemporai = — -r~ ' T ' cosh ^-1/1 + cos^ -. 

ncAn fsmo ^ 
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The coupling function k{z) has units of inverse length, therefore, the dimen- 
sional equivalent of ( |4.2CI| ) scaled by (|A.7a|) is 



nAn , 
A* = - — 1 - 



\UJ\ 



The defect width ( ^1211) scaled by is 



FWHM 



nAn 



k 



tanh 



-1 



'2\uj\Vuj^ + n^k^ + n^k^ - 
2nk 



With this information, we can construct the dimensional parameters de- 
scribing the numerical experiments in Section |^. Table |1] describes the defects 
and Table ^ describes the gap solitons. For simplicity, all pulse measurements 
are computed for f = so the spatial width must be used. 



Experiment 




k 


n 


A, (m-i) 


Defect FWHM (mm) 


1 


±1 


4 


1 


680 


1.6 


2.1 


-1 


2 


2 


680 


3.1 


2.2a 


-1 


1.6 


2.5 


680 


3.9 


2.2b 


-1 


1.33 


3 


680 


4.7 



Table 1: Experimental parameters describing the defects in numerical simu- 
lations. 



Experiment 


S 


/ (GW ■ cm-2) 


Gap Soliton FWHM (mm) 


1.1 


0.9 


490 


2.3 


1.2 


2 


1800 


1.3 


1.3 


n/2 


1300 


1.5 


2.1a 


0.9 


490 


2.3 


2.1b 


0.6 


230 


3.4 


2.2 


0.45 


130 


4.4 



Table 2: Experimental parameters describing the gap solitons in numerical 
simulations. 
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B Defect gratings with a prescribed mode 

In this section, we define a simple procedure for generating grating profiles 
with a given band gap, eigenvalue, and an eigenmode of prescribed shape. 
We may specify uj and look for solutions of the form (see (f4.4| )) 



g-ia;Tgja30(Z) 



f{Z) 



(B.l) 



where v± are constants, dzQ{Z) = V{Z), and f{Z) is a real scalar function 
such that f{Z) ~ e"''!^! as \Z\ oo. 



Then 



dE+ 



-iujT I el 



if' + iVf) e'^^^^v+. 



rt 

Letting g = j-i this may be rewritten 



dZ 



{g±iV)E+. 



Similarly Then ( [4.1| ) becomes 

{uj ± ig)E± + kE^ = 0; 

By m), 



Cv 



[iv + ig)e 



ie 



KC 



-ie 



ne 



ie 



{ijj — ig)e '® 
If (^) defines a solution of (^?T|) then 



det C = LO^ + g^ - = 0. 



Therefore, 

k{Z) = ^u^ + g^{Z). 
Note also that if f{Z) ~ e^'^''^' then \g\ ±kasZ =poo and 



lim k(Z) = Vuj^ + kP'- 



Z— »oo 



(B.2) 
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The width of the gap is then equal to 2koo- 
If w is a null eigenvector, then 



^ = -e-e^t±!^. (B.3) 



By (p.2|), we have that 



uj + ig 



1 



and therefore 



— = e*°, a real, 
f _ 



Examination of the eventual solution shows that a merely reflects the invari- 
ance of the equations under a constant phase shift. We therefore set a = 
in the remainder of the argument. Then, in order to satisfy ( p.3| ). 



on ~^ + ^3 
I'd = arg . 



K 



Since k is positive 



Then 



e = --arctan-. (B.4) 



We may use this method to construct defects that support defect modes 
of arbitrary shape with prescribed exponential decay. If we choose a function 
with different exponential decay rates as Z — > ±oo, then k±oo will take two 
different values. 

In the case a; = 0, we find the same discontinuous limiting behavior as 
in (|4.14| ). Fortunately, this is merely due to the inadequacy of the polar de- 



composition implicit in the definition of Q{Z) in ( [B.l| ). For = 0, the entire 
calculation may be repeated with = 0, and a smooth solution generaliz- 
ing (|4.16|) is generated: 

k{Z) = ±g{Z), (B.6) 
v_ = :jpiv+. (B.7) 
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Example 



The defect gratings of Section ^]3| which generahze the dark sohton gratings 
can be obtained as follows. Let 

/ = sech"(fcZ). 

Then 



k{Z) = y cj2 + tanr (kZ); (B.8) 

1 nkta.iah(kZ) _ , 
e{Z) = - arctan ^ ^; (B.9) 

2 LU 

V{Z) = ^^k^^ech^kZ) 

^ ' 2{uj^ + n^k^i8Ji\i\kZ)) ^ ' 



And 



, I i , nk tanh (kZ) 



Setting n = 1, we recover the "dark soliton defect" of Section [4.2| . By varying 
n, while keeping the quantity cu^+n^/c^ fixed, we generate a family of gratings 
of variable widths with identical band gaps. 

Numerical computations indicate that for n > 1, the system supports 
multiple eigenmodes obeying the following rule. For n > 0, the defect sup- 
ports a total of 2 [n] — 1 eigenmodes where \n \ is the smallest integer greater 
than or equal to n. The ground state has frequency Uq = u and spatial decay 
rate nk, whereas the excited states occur in pairs with spatial decay rates 
given by {n — j)k and frequencies: 

uj±j = ±^/uj^ + {2nj - p)k^ (B.12) 

for all 1 < j < n. It should be possible to derive this formula exactly 
using methods of complex analysis developed to study bound states of the 
Schrodinger equation with potential [Q, as well as expressions for the asso- 
ciated bound states. 
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